QUASISTATIC RHEOLOGY AND THE ORIGINS OF STRAIN. 
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RESUME : On confronte l'approche microscopique par simulation numerique discrete des materiaux 
granulaires de type solide a leurs proprietes rheologiques macroscopiques. On cite des systemes modeles 
dont les reponses, en deformation, un increment de contrainte different qualitativement, bien que la 
repartition des efforts Pequilibre soit tres similaire. Des resultats sur la sensibilite aux perturbations des 
reseaux de contact elastoplastiques permettent de distinguer deux regimes rheologiques, selon que leurs 
intervalles de stabilite, en termes de contraintes, se reduisent ou non zero dans la limite thermodynamiquc 
('fragilite' macroscopique) . On en evoque de possibles consequences. 

ABSTRACT: Features of rheological laws applied to solid-like granular materials are recalled and 
confronted to microscopic approaches via discrete numerical simulations. We give examples of model 
systems with very similar equilibrium stress transport properties - the much-studied force chains and force 
distribution - but qualitatively different strain responses to stress increments. Results on the stability 
of elastoplastic contact networks lead to the definition of two different rheological regimes, according to 
whether a macroscopic fragility property (propensity to rearrange under arbitrary small stress increments 
in the thermodynamic limit) applies. Possible consequences are discussed. 
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1 Scope 

This is a brief introduction to the rheology of solid-like granular materials in the quasistatic 
regime, with a special emphasis on the microscopic origins of strain, and on discrete numer- 
ical simulations of model systems. Rather sophisticated macroscopic phenomenological laws 
have been proposed [H [21 [3], but, in spite of many microscopic studies [H [5], with numerical 
tools [7] in particular, their relation to grain-level physical phenomena is not fully understood. 
Consequently, we mostly address basic, qualitative aspects on model systems. Moreover, we 
specialize on cohesionless, nearly rigid grains, and to small or moderate strain levels (excluding 
continuous, unbounded plastic flow). Despite the many insufficiencies of present-day modelling 
attempts, interesting directions for future research, elaborating on preliminary results, can be 
suggested. We recall a few basic concepts (section 2), some of the macroscopic phenomenology 
of solid-state granular mechanics (part 3), and the necessary elements of a microscopic model 
(section 4). Then, properties of simple model systems studied by numerical means, in the large 
system limit, are discussed both in frictionless (section 5) and in frictional (section 6) systems. 
Section 7 suggests broader perspectives and speculations. 



2 The constitutive law approach: basic ideas. 

On setting out to identify a constitutive law for a solid material, one has to rely on some 
postulates that are worth recalling in the context of granular materials. Such a law should 
locally relate stresses a to strains e , or, more appropriately for granular systems, stress rates 
& and strain rates e should determine each other for a given internal state of the system. This 
state is to be conceived of as specified once the values of some state variables (a finite number p 
of quantities, including the stress tensor itself, that exhaust the macroscopic description of the 
system) are known. One may write, at each time t: 

£(x,i) =T (i(x,i),£(x,t),{a(x,i)}) or |(x, t) = Q (&.(*, t),g(x, t), {a(x, i)}) (1) 

x standing for any 'point', in the sense of continuum mechanics, in the sample, i.e. a represen- 
tative volume element from the microscopic point of view, {a} is the set of unspecified state 
variables (ai)i<i< p . Their evolution should be ruled by similar equations: 

di(x,t) =w(|(x,t),ff(x,t),{a(x,t)}) . (2) 

Once eqns. Q] and [2] are gi ven > an d supplemented with the appropriate boundary conditions, one 
also needs the initial values of a and a to be able to predict the evolution of the system for, 
say, a prescribed history of stress. The prediction of the initial state of the system is in general 
beyond the scope of the rheological laws we are dealing with, as it is the result of a process 
that might involve rapid flow. (One exception is the construction of a sample under gravity by 
successive deposition of thin layers at the free surface. If the initial state of a freshly deposited 
layer is known, one may apply solid-like rheology to the rest of the sample, which deforms very 
little under the weight of the new layer, and thus, iteratively solving the appropriate boundary 
value problem, calculate the initial state of a whole system as a result of its construction history. 
This procedure can be applied to silos [8] and granular piles [9]). 

The suggestions, put forward in the recent literature [10] , to look for direct relationships 
between stress components that result from the construction history of a sample, are attempts 
to model the assembling procedure, rather than the response to stress increments. The proposed 
relations are not constitutive laws in the sense of eqns. [1] and [2] they ignore strains, and they 
are not local (they depend on sample shape and boundary conditions). 

One thus needs some a-priori knowledge of the initial stresses: however small the com- 
ponents of the stress tensor, the orientation of its principal axes and its level of anisotropy are 
important. Cohesionless grains do not spontaneously assemble in any 'natural state', once sub- 
mitted to some externally imposed stresses they form packings the structure of which depends 
on those stresses. Functions J 7 , Q and H of relations CD and [2] must be discontinuous at a = 0. 
In rheometric experiments one needs in principle to check for sample homogeneity. 

3 Macroscopic aspects 

Constitutive laws like eqns. [T] and [2] are studied in soil mechanics [3, 5, 6j. In order to extract some 
information on such laws from experiments, it is convenient to choose configurations in which 
stresses and strains are expected to be homogeneous. This leads to the design of rheometers, 
the most often employed one in soil mechanics being the triaxial apparatus, sketched on fig. [TJ 
Samples are submitted to axisymmetric states of stress, the axial stress a\, or the axial strain e%, 
are controlled via the relative motion of the end platens, while the lateral pressure p = a% = f 3 is 
exerted through a flexible membrane by a fluid. With some care (e.g., measuring strains directly 
on the sample in the central part away from the rigid platens) it is possible, with the most 
sophisticated devices, to record strains with an accuracy of the order of 10" 6 PJ H2]- Other 




Figure 1: Sketch of a triaxial experiment. 
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Figure 2: Schematic variations of deviator q (continuous curves) and volumetric strain e v (dotted lines) 
for a dense and a loose sand. The curve marked with arrows is observed on decreasing g or ei. 

rheometers [13] are the 'true triaxial' apparatus, which may impose three different principal 
stresses to a cubic sample, and the 'hollow cylinder' apparatus, which allows for rotation of the 
principal axes of stress and strain. In a typical triaxial experiment, one starts from a given 
state with, e.g., hydrostatic stress (a_ = pi). Then, most often, ei is increased at a constant 
(slow) rate, while lateral pressure p is maintained constant. Axial stress a± - or, equivalently, 
deviator q = a\—p- and the lateral strain 62 = £3 (or, equivalently, the relative volume increase, 
e v = — tr(e)) are measured. Evolutions of q and e v as ei monotonically increases are schematically 
represented on fig.[2j While the density and the deviator steadily increase, in a loose sample, until 
asymptotic constant values are reached, dense samples are initially contractant, then dilatant 
and the deviator curve passes through a maximum. From the beginning, the increase of q with 
ei is not reversible: if the direction of deformation is reversed, the same curve is not retraced 
back: the decrase of the deviator is steeper, with a slope comparable to that of the tangent at 
the origin of coordinates on fig[2j 

As €1 increase, the curves approach a plateau, corresponding to a final attractor that 
is called the 'critical state' [3] and deemed independent of the initial conditions (density and 
deviator should coincide for the loose and the dense sample on fig. [2] at large ei). However, the 
approach to this critical state is often hidden, in dense samples at least, by instabilities leading 
to strain localization in 'shear bands' (whose thickness is of the order of a few grain diameters). 
The development of these localizations was observed by X-ray tomography [14J . It is sensitive to 



sample shape and boundary conditions, but usually occurs in the vicinity of the observed 'stress 
peak'. Localization in loose samples, if it exists, is far less conspicuous. As sample homogeneity 
is lost, localization precludes the interpretation of rheometric tests in terms of constitutive laws, 
which should therefore be restricted to the 'pre-peak' part of the curve. 

An important feature of rheological tests on solid-like granular samples is their indepen- 
dence on physical time (see however, the remarks of section 7). Replacing t by <p(t) with a 
monotonically increasing function <p does not change stress-strain curves. Functions J 7 , Q , 7i 
should therefore be homogeneous of degree one in e or in &. Such tests are supposed to be quasi- 
static, as a sequence of equilibrium states is explored. Dynamical characteristics, such as masses, 
are regarded as irrelevant. Dry sands and water-saturated ones, provided static properties are 
the same, should exhibit the same behaviour. 

Eqns. [T]and [21 with due account for material symmetries, remain extremely general. The 
nature of internal variables a is not easy to guess if only stresses and strains are measured. The 
packing fraction is often chosen, because of its influence on the behaviour (fig. [2]). 

To account for irreversibilies and failure, the most commonly invoked laws are of the 
elastoplastic family. Those (although usually written differently) can be cast in the form of 
eqns. [T] and EJ We do not review elastoplasticity here. Its application to soil mechanics is 
presented in ref. [3]. Connections to limit analysis (calculation of a limit load beyond which 
unlimited plastic flow occurs) are discussed in [15]. Ref. [16] is a pedagogical introduction 
with examples of calculations with simplified laws. The complex behaviour sketched on fig. [2] 
requires quite sophisticated elastoplastic laws, with many parameters. Those laws should be 
non- associated |16| . and involve work hardening. Very roughly speaking, this means that the 
direction of plastic irreversible strains is not simply related to the failure condition on stresses, 
and that failure is a gradual process. A simplified law with correct qualitative properties in terms 
of global failure under monotonically varying loads [16] involves 4 parameters. A more elaborate 
and quantitative one, Nova's law [17] . requires 8 parameters. Incremental [1] or hypoplastic |18j 
laws are less traditional. They directly state relations like eqns. [1] and [2j J- and Q, although 
positively homogeneous of degree one, should be non-linear functions of e and & respectively, to 
account for irreversibility. This is directly postulated in such approaches, while elastoplastic laws 
describe such a behaviour through work hardening. Ref. [18] defines a simplified hypoplastic 
law with 5 parameters. 

Obviously, it is highly desirable to identify parameters with a physical meaning, connected 
with the microscopic mechanisms of deformation under stress. This would ease and guide the 
choice of a constitutive law, contribute to assess its range of validity (e.g., in terms of stress 
magnitudes), and reveal the influence of microscopic characteristics of a given material on its 
macroscopic behaviour. 

4 Ingredients of a microscopic model, discrete simulations. 

Numerical simulation methods are described in ref. [7]. They deal with simple models of granular 
materials, suitable to investigate the microscopic origins of constitutive laws. Here, we mainly 
discuss the simple cases of discs (2D) or spheres (3D). Three different kinds of ingredients are 
needed: geometric, static and dynamical ones. First, grain shape and polydispersity have to 
be specified. Then, static parameters are those defining equilibrium contact laws. In the most 
simple model, a normal stiffness constant Kn expresses a proportionality between normal force 
/at and normal deflection h of a contact (which is modelled as an interpenetration depth), a 
tangential stiffness Kt incrementally relates tangential contact forces fx to tangential relative 
displacements, and Coulomb's condition \/t\ < h/n should be satisfied, with a friction coefficient 
fx, sliding being allowed (such that the work of fx is negative) if it holds as an equality. One may 
regard stiffness constants as a mere computational trick to forbid grain interpenetration. They 



may also be chosen with correct order of magnitudes, comparing them to typical estimated values 
of dfisr/dh in contacts, under some given stress level. It can even be attempted to implement 
accurate contact laws, such as the Hertz-Mindlin-Deresiewicz [19] ones for smooth elastic spheres 
with friction. 

Finally dynamical parameters are related to inertia (masses, moments of inertia) and 
kinetic energy dissipation (e. g. viscous damping in contacts). In slow, quasi-static evolutions, 
those parameters should be irrelevant, the behaviour should be determined by the geometric 
data, along with /x, ratio Kt/Kn and parameter K^/P (in 2D) or K^/{Pd) (in 3D), which 
measures the deflections of contacts relative to grain diameter d under typical forces [23] . 

The most widely used simulation methods [7j, molecular dynamics (MD), or contact 
dynamics, rely on time integration of dynamical equations of motion. They have been used 
to simulate biaxial (in 2D) [20] or triaxial [21] tests. Such calculations are usually made at 
constant axial strain rate, assuming the system remains close enough to equilibrium at any 
time for the evolution to be regarded as quasi-static. They cope with a few hundreds or a 
few thousands of grains. They successfully produce stress-strain curves whose broad features, 
on the scale of e ~ 10%, are those of fig. [2j For instance Thornton's simulations [2lj yield 
a maximum deviator criterion that coincides with some experimental observations. However, 
stress-strain curves, for given loading histories, are still rather noisy on a smaller scale (e ~ 1%), 
especially as the peak is approached (strain values corresponding to the peak deviator are not 
accurate). It seems necessary to study the form of such curves and investigate the regression 
of fluctuations with greater care, for two reasons: first, one might wish to obtain accurate 
estimations of rheological laws in the macroscopic limit; then reliable numerical data on the 
effect of perturbations on equilibrium states would give insight on the microscopic mechanisms 
of deformation. Those should be accounted for in theoretical attempts to relate rheological laws 
to grain-level phenomena. The next sections are therefore devoted to biaxial compressions with 
2D systems of disks, in which deviator q = o\ — o~i (keeping the notations introduced for the 
triaxial test) is stepwise increased and one studies the effect of small stress increments imposed 
on equilibrium configurations. 

5 Response to stress increments: frictionless grains 

Assemblies of frictionless grains are particularly appealing because of two remarkable properties, 
that are established and discussed in ref. [22]. The first one, the absence of hyperstaticity in the 
limit of large contact stiffness, means that the contact network is barely sufficient to support 
stresses, and that the sole condition that only closed contacts can transmit a force, along with the 
force (and torque) balance equations, is sufficient to calculate all contact forces. The second (the 
standard mechanical energy minimization property) states, for rigid grains, that the potential 
energy of external forces has to be minimized under the constraints of no interpenetration. It 
allows to discuss the stability of equilibrium states. 

Together, both properties entail [22] that force-carrying structures in assemblies of rigid 
frictionless and cohesionless disks (in 2D) or spheres (in 3D) are isostatic: there is a one to 
one correspondence between external forces and contact forces, and between relative normal 
velocities in the contacts and grain velocities. We exploited this [23J to study the response of 
disordered systems of disks to stress increments, by a purely geometric procedure we called the 
'geometric quasi-static method' (GQSM). In such systems, the stress-strain curve is a staircase 
and the procedure tracks the elementary steps. In stability intervals, the rigid contact structure 
supports the stress without motion, until one contact force becomes negative. This contact has 
then to open, initiating a rearrangement, hence a strain increment. Motion stops when another 
contact closes and a new stable equilibrium is reached. In fact, this might require several contact 
replacements, which are operated one by one in this algorithm. On opening one contact, the 
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Figure 3: Velocity field at the beginning of a rearrangement in a system of 1024 disks in a square box 
undergoing biaxial compression. The contact between the disks marked in gray (bottom right) is opening. 



ensuing velocity field is, up to a positive factor, geometrically determined. One example is 
displayed on fig. [3l Those fields form complex vortex patterns extending through the whole 
sample. Displacement fields corresponding to small strain increments have the same aspect. 

Statistics of stress and strain steps for the beginning of a biaxial compression (close to 
the isotropic state of stress) were studied [231 HI] • It was found that intervals of stability 8q 
are exponentially distributed and scale with the number of disks N as N~ a with a ~ 1.1, 
while 'axial' (conjugate to the largest principal stress that is being incremented) strain steps 
are power-law distributed, the density function decrasing as (<5e) _ ^ 1+M \ with \i ~ 0.5 for large 
values, and scale as with j3 ~ 2.1. 

As stability regions in stress space dwindle to nothing in the thermodynamic limit (a 
macroscopic 'fragility' property [251 I22j). equilibrium states will be rather elusive: it is im- 
possible, in a real experiment, to control stress levels with perfect accuracy. Although each 
equilibrium configuration is rigid, any level of noise in a macroscopic system should generate 
fluctuations, because rigid configurations become unstable, and the system should keep visiting 
several equilibrium states. 

Since the power law distribution of strain increments does not admit a mean value, the 
accumulation of successive 5e steps generates a Levy process [26], and the strain step corre- 
sponding to a given deviator increment remains impredictable. Although, as N increases, the 
typical size of steps decreases, the staircase stress-strain curve does not become smooth because 
of the statistical importance of large strain increments. It was also observed [23] that the statis- 
tics of strain variations corresponding to given stress intervals do not depend on N. Moreover, 
this distribution is the same for the GQSM and for a more conventional MD calculation (with 
nearly rigid grains, K^/P = 10 5 ). As MD results, introducing additional static and dynami- 
cal parameters, are statistically indistinguishable from GQSM ones, one may conclude that the 
mechanical response is determined by geometry alone. These results are summarized on fig. 0] 

This very singular behaviour - a constitutive law cannot be defined - calls for additional 
investigation of its origins and range of validity. It is worth pointing out that the observation 
of the large 5e values and the Levy-stable distribution are due to the rearrangements in which 
several (sometimes many) contacts have to be replaced. If we now adopt the approximation 
of small displacements (ASD) in which displacements away from a reference configuration are 
dealt with as infinitesimal, and hence normal unit vectors between neighbouring disks are kept 
constant, then it can be shown [22] that only one contact replacement will be enough to reach 
the next equilibrium state. Then, the distribution of strain steps admits a first and a second 
moment, and the staircase, within the ASD, should approach a smooth curve. 
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Figure 4: Distribution of axial strain increments corresponding to fixed Aq = 10~ 3 p measured for 4 
sample sizes (up to N = 4900) with GQSM and 2 sizes (1000 and 3000) with MD. Note the indepen- 
dence on the numerical method (for large values) and the absence of convergence, as N increases, to a 
deterministic response. 

In fact, the ASD with frictionless grains has deep consequences: it entails |22] the unique- 
ness of equilibrium states. There is no need to obtain a stress-strain curve in an incremental way, 
as stresses and strains are in one-to-one correspondence. Uniqueness also implies the absence 
of irreversibility, since there is no history dependence. Assemblies of slightly polydisperse disks 
placed on a regular triangular network were studied within the ASD (with the unperturbed 
lattice as reference). The approximation is well controlled in that case because of the small 
polydispersity parameter. This system might be called elastic (although disks are rigid). Its 
behaviour was shown [27\ [22] to be analogous to that of a mobile point requested to stay within 
a convex part T>, limited by a smooth surface S, of a three-dimensional space (the analog of 
the set of values permitted by impenetrability conditions in strain space). Once submitted to 
an external force (its 3 coordinates are the analogs of the stress components) its equilibrium 
position is the point of £ where the tangent plane is orthogonal to the force. Hence a smooth 
correspondence between force and displacements. Upon incrementing the force, the displace- 
ment is inversely proportional to the curvature of E. As a consequence of these properties [22] . 
the macroscopic response (displacements) to some small localized force superimposed on a pre- 
imposed stress field (Green's function) is, for this model, the solution to an elliptic boundary 
value problem (akin to elastic problems for incompressible materials). Within the ASD, finding 
an equilibrium state amounts to solving a convex minimization problem. This is also the case, 
without the ASD, for networks of cables, for which the same kind of elasticity applies. 

Let us summarize the main conclusions of these studies on rigid frictionless grains. 

(i) The stress-strain curve is a staircase with phases of stability, with just enough contacts 
to carry the forces, alternating with rearrangements, with just enough contact openings to allow 
some deformation. Those rearrangements are non-local events, involving the whole system, (ii) 
Any macroscopic stress perturbation causes some rearrangement. (Hi) No deterministic con- 
stitutive law applies to disordered assemblies of rigid, frictionless disks, (iv) The mechanical 
response is determined by the sole geometric data, (v) Different systems might exhibit equilib- 
rium states with extremely similar properties (in terms of force distributions) and both satisfy 
properties i and ii. However their mechanical behaviour might be drastically different: conclu- 
sion Hi applies to disks without the ASD, whereas disks within the ASD or cable networks abide 
by some form of elasticity, rearrangements being reversible. 




6 Response to stress increments: grains with friction 

Let us now report some results on disordered systems of nearly rigid (K^/P = 10 5 ) systems 
of disks with a friction coefficient /i > in the contacts [28, 24J. Just like in the frictionless 
case, one may either try to track elementary stability intervals and rearrangements, or resort 
to molecular dynamics (introducing additional parameters). However, the specific properties 
of frictionless systems are lost: the contacts may transmit tangential forces and the network is 
hyperstatic; potential energy is no longer minimized at equilibrium. To discuss the stability of a 
given contact network, one needs to introduce elasticity in the contacts. One can then perform 
static elastoplastic calculations: the system is to be regarded as a network of springs, plastic 
sliders and no-tension joints, the displacements and rotations of the disks can be computed 
for each applied stress increment, via an iterative process [29]. Such studies of elastoplastic 
networks [291 [30] with static methods are surprisingly rare, especially in comparison to the vast 
literature on numerical simulations of elastic networks and brittle fracture (see e.g. |31j) The 
systems studied in [28] were prepared without friction, and are thus very dense. Several simple 
sizes were studied (N ranging from 1000 to 5000), and it was found that the initial contact 
network, corresponding to isotropic stresses at the beginning of the biaxial compression, is able 
to support a considerable stress deviator (q/p = 0.81 ± 0.06 for \i = 0.25) in the large system 
limit. Therefore, such dense systems with friction are not fragile in the sense of section 5. 
Stress-strain curves for /x = 0.25, Kt/Kn = 0.5 for the beginning of the biaxial compression 
are shown on fig. [5j In this regime, that we call strictly quasi-static, the curve is smooth, 
the successive equilibrium configurations form a continuum. The scale of strains is set by the 
stiffness constants in the contacts. The behaviour is inelastic and irreversible from the beginning 
of the biaxial compression, as the proportion of sliding or opening contacts steadily increases. 
Eventually, some instability occurs, the initially present contacts can no longer support the 
stresses. Interestingly, this appears to happen for a deviator value that does not sensitively 
depend on stiffness constants [24] . To proceed further (as the current state of the static algorithm 
does not clearly determine one direction of instability and no analog of the GQSM is available) , 
we resorted to molecular dynamics. Successive equilibria corresponding to stepwise increasing 
deviator values were obtained, and a staircase-like stress-strain curve was observed, signalling 
the frequent occurrence of instabilities and strain jumps (fig. [5]). To check whether a well- 
defined stress-strain relation is approached as N — > oo, averages and mean standard deviations 
of q(e\) for < e\ < 0.02 were computed with several samples of different sizes, and the results 
(fig. [6]) do indicate that a smooth curve is approached. Similar results are obtained for e v versus 
e\. Any strictly quasi-static interval appears as a vertical segment on fig. [6j The existence 
of a limit for large N requires in fact the fragility property to apply in the staircase regime. 
However, it is expected that although any positive increment 5q will entail a rearrangement in 
a macroscopic system, contact structures will withstand finite negative Sq's. Fig. [5] shows that 
relatively large q intervals could be accessed by static calculations, from intermediate equilibrium 
configurations in the staircase regime, upon decreasing q. Moreover, it is observed that many 
contacts stop sliding on reversing the motion. It would be interesting to investigate the response 
to differently oriented stress increments, and to delineate the (history-dependent) strictly quasi- 
static, non-fragile domain around a given equilibrium state. Taking the mobilization of friction 
into account - i.e., replacing Coulomb's inequality by an equality for all sliding contacts - 
it can be observed that the indeterminacy of forces is greatly reduced within the staircase 
regime in the monotonic biaxial compression [201 [28] , which suggests an analogy with isostatic 
frictionless systems. Grain motions, rearrangements and spatial distribution of strains were 
studied by Williams and Rege [32] , and by Kuhn [33] . Similar patterns as those of fig. [3] were 
observed. Thin non-persistent (unlike shear bands) 'microbands' concentrating the strain were 
also reported [33]. 






Figure 5: Main plot: aspect of deviator/ axial strain curves for e\ < 0.02 (MD calculations). The 
dotted lines are results obtained with the static method on going backwards on the same stress path (q 
decreases). The inset - note the blown- up e scales - shows q and e v versus ei in the strictly quasi-static 
regime 
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Figure 6: q{e\) for 3 sample sizes: the larger N, the darker is the shaded zone extending to one standard 
deviation above and below the average. The inset displays the average standard deviation throughout 
the ei range, for q and for the volumetric strain, versus N~ x / 2 . 
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Figure 7: Axial stress vs. axial strain curve, from ref. 0D]. The experiment was stopped at constant 
stress several times, hence the creep intervals. Note the stiff response when the constant rate compression 
is resumed. 

There is therefore some evidence that conclusions (i), (ii), and (iv) of section 5 are still 
valid in systems with friction within the fragile 'staircase ' regime, the essential differences being 
the role of the friction coefficient itself (the behaviour appears to be essentially determined by 
the geometry and fj), the existence of a well-defined macroscopic stress-strain curve and that of 
strictly quasi-static regimes within which a given network of contacts is able to support a finite 
stress range, all strains being due to the finite stiffness of the grains themselves. 

7 Some perspectives and speculations. 

The non-local aspect of rearranging events, reflecting the strong steric hindrance in dense pack- 
ings of impenetrable bodies, could be expected to preclude the definition of a local law. How- 
ever, one might think that a great number of such long-distance correlated motions of very 
small amplitudes could, once aggregated, build up a strain field devoid of long-range corre- 
lations. Specifically, strains could be localized on some 'microband' pattern (as reported by 
Kuhn |33j ) during one elementary rearrangement, but the random superposition of lots of such 
non-persistent structures could destroy the long-distance correlations. Similar ideas were fol- 
lowed by Torok et al. [34] : their schematic model assumes macroscopic deformation to result 
from an accumulation of slides on temporary slipping surfaces. It would be interesting to test 
such a scenario, which requires an accurate numerical computation of directions of instabilities 
of granular assemblies with a given contact network. That 'deformation consists in a number 
of arrested slides' is an idea already put forth by Rowe in 1962 (ref. [35j,p. 514). His classical 
'stress-dilatancy relation' [35J could thus be founded on a microscopic analysis. 

In parallel to the analysis of stress-strains relationships we have been reporting here (we 
focussed on eqn. [TJ, microscopic studies have tried to define internal state variables of granular 
systems and to relate them to stresses and strains (i.e. suggesting a form of eqn.[2|). For instance, 
the density of contacts and some parametrization of the distribution of their orientation (called 
fabric or texture) have been studied, their evolution can be related to strains [36], [37] (eqn. [2]) 
and their values can be correlated to the possibly supported stress orientations |38j (role of a in 
eqn. [Tj. A brief presentation of the possible use of packing fraction and fabric as work-hardening 
variables in a plasticity theory is ref. [39|. 




Finally, let us briefly speculate about possible consequences of the existence of strictly 
quasi-static and fragile regimes. Sand specimens were observed to creep: under constant 
stresses |40[ |4"T] (e.g., on stopping a triaxial test and maintaining constant stresses) strains 
vary very slowly, over hours. (Let us quote ref. [35] again: 'The time to equilibrium increases to 
many days as the peak strength is approached'). A new equilibrium might be approached, which 




can be rather distant. When the slow controlled strain rate is resumed, the response to both 
positive and negative q increments is quite stiff. With due account to possible aging phenomena 
in the contacts [12], it is tempting to suggest the following explanation: once left to wait at 
constant stress, the grain pack, which is highly sensitive to noise, slowly drifts in configuration 
space, until it reaches a state with a finite stability range (hence a stiff response of the 'strictly 
quasi-static' type). There appears thus to be interesting connections between slow dynamics 
and fragility. 
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